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ABSTRACT 

New images of M31 at 24, 70, and 160 /xm taken with the Multiband Imaging Photometer for Spitzer 
(MIPS) reveal the morphology of the dust in this galaxy. This morphology is well represented by a 
composite of two logarithmic spiral arms and a circular ring (radius ~ 10 kpc) of star formation offset 
from the nucleus. The two spiral arms appear to start at the ends of a bar in the nuclear region 
and extend beyond the star forming ring. As has been found in previous work, the spiral arms are 
not continuous but composed of spiral segments. The star forming ring is very circular except for a 
region near M32 where it splits. The lack of well defined spiral arms and the prominence of the nearly 
circular ring suggests that M31 has been distorted by interactions with its satellite galaxies. Using 
new dynamical simulations of M31 interacting with M32 and NGC 205 we find that, qualitatively, 
such interactions can produce an offset, split ring like that seen in the MIPS images. 
Subject headings: galaxies: individual (M31, M32, NGC 205) — galaxies: spiral — galaxies: ISM 



1. INTRODUCTION 

M31 is the nearest (d ~ 780 kpc, iStanek &: Garnavicbl 
1998 [lRich et all2005fl external spiral galaxy fSb. lHubblel 
1929) making it one of the premier laboratories for un- 
derstanding spiral galaxies like our own. Tracing the spi- 
ral structure in M31 is an obvious first step, which has 
been done with varying degrees of success for many years. 
However, the high inclination (74-78°) of M31 makes it 
difficult to determine its large-scale structure. In particu- 
lar, there is no clear consensus regarding the organization 
of spiral structure in this important galaxy. The earliest 
quantitative studies focused on optical tr acers such as 
HII re gions l| Arp||1964|:|Baa4e Arp||1964fl. OB associa- 
tions (lHodgelll979l : Ivan den Bergblll991j) . and dark neb- 
ulae ijHodgdl 1980(1 . These studies generally found struc- 
tures consistent with a two-armed spiral (f or the case for 
a one-armed spiral see lSimien et all 119781) . but invoked 
disturbances such as varying inclinations and warps to 
account for the deviati ons from simple logarithmic spi- 
rals. More recently HI llGuiberdll974HBraunlll991^ and 
CO l|Loinard et all 119991: iNieten et all 12005ft have been 
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used. The velocity information available for these gas 
tracers can be used to disentangle the pileup of informa- 
tion due to M31's high inclination. These studies have 
found strong evidence for a two-armed spiral structure, 
but still suffer from the difficulty in connecting spiral arm 
segments into a coherent pattern. 

The difficulty in finding coherent spiral structure in 
M31 is plausibly related to the interactions of M31 with 
its two nearest satellites, M32 and NGC 205. Dis- 
turbances in M31 are revealed by tidal features far 
outside M31's disk visible in de ep, wide-field optical 
imaging (llbata et all l200ll 120051: iFerguson etaD 120021: 
iLewis et alJ I2004|L M32 is understood to be the per- 
turbcr of the spiral arms llBvrdlll978l Il98l while NGC 
205 has been modeled as th e cause of the warp s een in the 
optical and HI disk of M31 <|Sato & Sawall986[) . Quanti- 
fying the effects of M31's satellite galaxies is complicated 
by the lack of measured p roper motions and uncerta inties 
in relative distances (e.g. iMcConnachie et al.l f2005K 

Another difficulty in quantifying M31's spiral structure 
has been the lack of appropriate tracers in the innermost 
regions of the galaxy. Optical tracers are overwhelmed 
by the bright stellar bulge, while HI and CO, which 
woul d be easily seen through the bulge, are larg ely ab- 
sent ijBrinks fc Shanejll984HMelchior et al.ll2000l i as the 
interstellar medium (ISM) i n the central region is com- 
posed mainly of ionized gas l)Devereux et al.lll994D . The 
far-infrared (far-IR) emission can circumvent these diffi- 
culties as it traces the dust in M31, which penetrates the 
stellar bulge and should trace all phases of the ISM. 

The first far-IR images of M31 were obtained with the 
Infrared Astronomical Satellite (IRAS) at 12, 25, 60, and 
100 /im. The most prominent feature in these images 
is the bright, ~10 kpc radius, ring of star formation 
llHabing et all Il984t iRicd Il993[) seen pr eviously in Ha 
images ljArplll964l iDevereux et al.lll994D . The Infrared 
Space Observatory (ISO) imaged M31 at 175 /im finding 
both the bright r ing and a fainter ring at larger radii 
l|Haas et al]ll998fl . The Midcourse Space Experiment 
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Fig. 1. — The M31 MIPS images are shown. The images are scaled using a quadratic root (a; 25 ) between -0.1 — 3 MJy/sr (24 fira), 
0.1 — 35 MJy/sr (70 ^m), and 0-90 MJy/sr (160 (im). The bottom image gives the 3 color composite of the MIPS images (with scaling 
adjusted to emphasize the differences between bands); the color gives an indication of the dust temperature. The images are oriented with 
the nominal position angle of M31 (38°) horizontal and are cropped to a common size. The nearly- vertical streaks present in the 70 and 
160 (im images are residual Ga:Ge detector instrumental signatures. M32 is clearly detected at 24 /im, but not in the other two MIPS 
bands. NGC 205 is detected in all three MIPS bands, but is located NW of M31 beyond the boundary of the cropped images. The two 
spots near the nucleus marked on the 160 fim image are discussed in the text. 
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Fig. 2.— The MIPS global fluxes of M31 are p lotted 
along w ith previous measurem ents f rom IR AS I Rice ct al .1 119881) . 
COBE IQdenwald et ai]ll99l) . ISO IHaas et alj|19981) . and MSX 
(Kracmcr ct al. 2002). A simple stellar population plus dust grain 
model (Marleau et al., in prep.) has been fit to the data giving 
a total IR luminosity of 1.7 X 10 43 ergs s — 1 which corresponds 
to a star formation rate of 0.75 Mq yr -1 using the calibration of 
IKennicuttl I199SI) . 



(MSX) provided an image of M31 at 8.3 /im revealing 
possible spiral structure inside the ~10 kpc ring. The 
resolutions of the IRAS (~ 1') and ISO (lf3) images re- 
veal the overall morphology, but are not high enough to 
reveal the structure of the spiral arms. 

We present new infrared images of M31 taken with 
the Multiband Im aging Photometer for Spitzer (MIPS, 
IRieke et al.ll2004|) instrument on th e Spitzer Space Tele- 
scope (Spitzer. iWerner et alJfeOO^ . The MIPS spatial 
resolution at the distance of M31 is 6"/23 pc, 18"/68 pc, 
and 40"/151 pc at 24, 70, and 160 /im, respectively. 
This paper concentrates on the newly-revealed details of 
M31's infrared morphology as shown in the MIPS images 

2. DATA 

The MIPS images of M31 were taken on 25 Aug 2004. 
A region approximately 1° x 3° oriented along the major 
axis of M31 was covered with 7 scan maps. Each map 
consisted of medium rate scan legs and cross scan offsets 
of 148" with lengths varying between 0?75 and 1?25 to 
efficiently map M31 a nd its satellite galax y NGC 205. 
The MIPS DAT v2.90 llOordon et al.ll2005l) was used to 
process and mosaic the individual images. At 24 /jm ex- 
tra steps were carried out to improve the images includ- 
ing readout offset correction, array averaged background 
subtraction (using a low order polynomial fit to each leg, 
with the region including M31 excluded from this fit), 
and exclusion of the 1st five images in each scan leg due 
to boost frame transients. At 70 and 160 /im, the extra 
processing step was a pixel dependent background sub- 
traction for each map (technique same as 24 /im). The 
background subtraction does not remove real M31 emis- 
sion as the scan legs are nearly parallel to the minor axis 
with the background regions being far from M31's disk. 
The final mosaics have exposure times of ^100, ^40, and 
~9 seconds/pixel for 24, 70, and 160 /im, respectively, 
and are shown in Figure ^ 

The IR spectral energy distribution of M31 is shown in 
Fig.0 The global MIPS fluxes measured in a 2?75 x 0?75 



rectangular aperture (the background was subtracted 
in the data reduction) are 107 ± 10.7, 940 ± 188, and 
7900 ± 1580 Jy for 24, 70, and 160 £im, respectively. The 
uncertainties are dominated by the systematic uncertain- 
ties in the MIPS calibration l|Rieke et alJl2004j) . There is 
excellent agreement between the MIPS and previous IR 
measurements . 

3. MORPHOLOGY 

The MIPS images in Fig. ^ are dominated by the 
~10 kpc ring. The deprojected ISO image showed that 
the ring is fairly circular, but with an obvi ous splitting 
near the position of M32 ijHaas et al.lll998|) . This split- 
ting is now seen to extend over ~l/4 of the ring. The 
outer ring seen with ISO is seen in the MIPS images, but 
appears to be relatively incomplete. Overall, the mor- 
phology of the disk of M31 in the infrared appears more 
disturbed on the left-hand (southwest) than right-hand 
(northeast) side of the galaxy. 

The MIPS images show emission extending beyond 
their edges along the major axis, implying that the in- 
frared extent of M31 is larger than 3°/40.9 kpc. A com- 
parison of the MIPS 160 /im image with the latest HI im- 
age (Braun et al., in prep.) reveals dust emission where 
there is HI emission. The gas and dust seem well mixed 
even in the outer regions of M31. 

The infrared emission in the nuclear region has a com- 
ponent with hotter dust than the rest of M31 (but 
also a cold component due to the emission seen in the 
160 /im image). This is consistent with the nuclear re- 
gion gas being domina ted by ionized gas as seen in Ha 
ijDevereux et aH ll994|) as well as emission from AGB 
stars in the bulge. As the flux from the nucleus fades 
towards longer wavelengths, the two spots above and be- 
low along the minor axis gain prominence. The p romi- 
nence of the spot above the nucleus was no ticed bv lRicd 
( 1993) and confirmed bv lHaas et alJ l)1998j) . These spots 
appear to be the ends of a bar and the beginning of spi- 
ral arms. The location of these spots and appearance of 
the nuclear region (especiall y at 70 /xm) is s imilar to the 
triaxial bulge simulations of Berman] <)2001jh 

To investigate the morphology in more detail, the im- 
ages have been deprojected assuming an inclination of 
75°, and displayed in Fig. EI Usually the inclination is 
picked visually to make the galaxy most spiral-like or 
circular in the deprojected frame. To be more objective, 
the inclination we use is the median value of the indepen- 
dently calculat ed inclinations of the spiral arm segments 
investigated bv lBraur] l)1991|h 

Arc-like structures dominate the morphology of the re- 
gion inside the ^10 kpc ring and appear to be spiral 
arms with starting points at ends of a bar. The exis- 
tence of a bar in M31 has been suggested in the past 
ijStark fc Binnevlll994|) . Since circularly symmetric ob- 
jects will deproject into elongated bars, the evidence for 
a bar is based on the prominence of the two spots near 
the nucleus at 160 /im, not on the appearance of a bar 
in the nuclear region at 24 fim. 

To determine the parameters of any spiral arms, Fig.0] 
gives plots in polar and Cartesian coordinates of the loca- 
tions of compact IR objects seen in all three MIPS bands. 
These regions are likely to be star forming regions (e.g., 
HII regions and OB associations) which are excellent spi- 
ral arm tracers. Point source photometry was done using 





Fig. 3. — The MIPS images deprojected using an inclination of 75° are shown in the left column. The right column shows the 24 fim 
image with the spiral arms shown in Fig. Hloverplotted. the M31-M32 interaction model (M32 as a blue dot), and M31 quiescent model 
(see gj. The M31-M32 interaction model image gives the true location of M32. If deprojected with the same assumptions as the images, 
the location of M32 would shift off the image (within the ~10' model insensitivity to the true position of M32). The model images give the 
integrated gas density which is nominally proportional to the dust given a constant gas-to-dust ratio (except in the inner region where it is 
clear that the gas-to-dust ratio likely varies significantly). The orientation of the images is the same as in Fig. and the size is 3° X 3°. 
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StarFinder (|Diolaiti et alJl2000j) after convolving the 24 
and 70 /zm images to the 160 ^m resolution using custom 
convolution kernels (Gordon et a!., in prep.) ensuring the 
same physical regions were measured in the three MIPS 
bands. 

In polar coordinates, a logarithmic spiral is a line with 
a non-zero slope and a circle is a horizontal line. The 
strongest feature is the wavy line around ln(radius) = 
3.75 representing the ~10 kpc ring. This wavy line is 
fit well by a circle offset from the nucleus by (5.5', 3.0') 
with a radius of 43'/9.8 kpc, except for the region around 
200° where the ring splits. The remarkable circularity of 
this ring makes it difficult to explain this structure as 
spiral arm segments. Two logarithmic spirals are plotted 
which roughly follow the spiral structure. These trailing 
spirals have pitch angles of 9° and 9? 5, and are phase 
shifted by 166°. These two spirals trace most of the 
structure inside the ring as well as some of the struc- 
ture outside the ring. These two spirals are not unique 
solutions, but reasonable given the disturbed nature of 
M31. The pitch angles used here are larger than found 
in previous work, consistent with not including the offset 
~10 kpc ring as part of the spiral structure. We found 
that our simple logarithmic spirals plus an offset circle 
give as g ood or better representation of the MIPS sources 
than the lBraunl (|1991|) spiral arm segments when the seg- 
ments were deprojected in the same manner as the MIPS 
sources. 

Finding that M31 has an offset, star forming ring which 
is not just a superposition of spiral arm segments raises 
the question: "What is the origin of the ring?" The ring 
could be due to internal resonances which can cause three 
differ ent types of rings: nuclear, inner, and outer iJButal 
1986). The M31 ring is clearly not nuclear or an inner 
ring which is directly associated with bar. The M31 ring 
could be an outer ring which is associated with the outer 
Lindblad resonance/radius of corota tion. The rad ius of 
corotation is placed around 18 kpc bv lBraunl (^991). The 
ring could also be caused by interactions like thos e seen 
in the Cartwheel galaxy ijHernauist fc WeTH fl993) . An 
interesting clue to the possible origin of the ring is the 
splitting of the ring by M32's passage through M31's disk 
(2J). Is this a coincidence or cause of the ring? A detailed 
study of the origin of M31's ring is clearly needed. 

While the logarithmic spirals plotted in Fig. 0] repre- 
sent the data well, there are clear deviations from these 
spirals in addition to the well defined circle of star for- 
mation. The known interactions between M31 and its 
satellite galaxies M32 and NGC 205 imply these inter- 
actions may be the cause of the deviations from regular 
spiral structures. 

4. DYNAMICAL MODELS 

The dynamical effects that M32 and NGC 205 
have on M31 were investigated by running a number 
of numerical simulations of the M31- M32 and M31- 
NGC 2 05 intera c tions u sing GADGET2 ijSpringel et all 
I200H ISpringell 12005]) . These simulations self- 
consistently model gravity and hydrodynamics. In ad- 
diti on, we have imp lemented star formation follow- 
ing IKatz et all (119961) . cooling using the MAPPINGS 
III (http://www.mso.anu.edu.au/~ralph/map.html) col- 
lisional ionization equilibrium models for [Fe/H|=-0.5, 
and a combination of thermal and kinetic feedback from 



supernovae (N avarro & White! Il993l with f v = 0.1) 
using the instantaneous recycling approximation. We 
used only the tree algorithm when calculating the grav- 
itational forces as we found that the default TreePM 
method introduced unacceptably large numerical arti- 
facts. 

The initial conditions for M31 w ere the "M31a" equi- 
librium disk-bulge-halo model of iWidrow fc Dubinskil 
( 2005), with 132092, 48000, and 357192 particles in each 
component respectively. A randomly-sampled fraction 
of the resulting disk particles were turned into gas par- 
ticles such that the radial gas density profile followed 
the sum of the HI (Braun et al., in prep.) and CO 
ijNieten et alJl2005|) gas profiles except in the inner re- 
gion where we assume a constant density to account for 
the mass of ionized gas. M32 and NGC 205 are modeled 
as poi nt particles, with radial velocities taken from lMateol 
( 1998). The current line-of-sight distances to the satellite 
galaxies were va ried within the unce r taintie s of current 
estimates (M32 : IMcConnachie et all l|2005l) : NGC 205: 
iFerrarese et al.l ((2000)), and the velocities in the plane 
of the sky were varied freely in order to find configura- 
tions where each satellite has recently passed through the 
diskofM31. 

In Figure [3] we show an image of one such simulation of 
M31 after an interaction with M32, assuming that M32 
currently lies 8 kpc behind the nucleus of M31, has a 
velocity in the plane of the sky of 200 km s _1 SSE, and 
has a mass of 1 x 10 10 M w (3-5 times larger than cur- 
rent estimates l|Mateolll99|) L We also show an image 
of a simulation with no satellite. The passage of M32 
through the disk of M31 results in a burst of star forma- 
tion that propagates outward through the disk resulting 
in a large hole similar in size and location to the observed 
splitting of the star forming ring. The location and age of 
the star formin g region NGC 206 is consistent with this 
star formation (jHunter et al.lll996|) . We find that such a 
feature rapidly distorts due to differential rotation, and 
so its present symmetric appearance indicates that the 
disk passage occurred very recently (20 Myr ago in this 
model). I f we use a mas s of 2 x 10 9 M Q for M32, as sug- 
gested bv iMateol (|1998f) . very little effect is seen on the 
disk of M31, suggesting that M32 may be more massive 
than previously thought. It is interesting to note that 
the input high density gas ring is pulled off-center (~1') 
by the satellite interaction, similar but smaller than the 
observed offset (6.3'). The simulated ring fragments into 
spiral arms and, therefore, is not as apparent as the ob- 
served ring. 

The mass of NGC 205 is uncertain {Mated i|1998D cal- 
culates a mass of approximately 8 x 10 8 Mq) and there- 
fore its ability to strongly affect the disk of M31 is equally 
uncertain. We have run a number of simulations of 
NGC 205, and while we can find orbits for which it passes 
through and temporarily distorts the spiral structure in 
the disk, we have been unable to find one which produces 
a hole that matches the observed split in the ring. How- 
ever, an exhaustive search of the parameter space has 
not yet been carried out. 

It is a generic feature in all of our models (with or 
without satellite interaction) that the gas in the central 
region of M31 ends up at a temperature greater than 
10 4 K, consistent with the observation of ionized gas in 
the nuclear region. This behavior is not seen in models 
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Fig. 4. — Polar and Cartesian plots showing the locations of point sources detected in all three MIPS bands at 160 (im resolution. The 
symbols have been scaled by the measured 160 /an fluxes. These plots are overlaid with two simple logarithmic spirals (green, one solid 
and one dashed) and an offset circle (black). 



with a fully exponential gas disk; therefore, the hot nu- 
cleus is a direct result of the central depression of the gas 
density. As central depressions are commonly seen i n the 
HI profiles of disk galaxies (e.g. iHewitt et alJll983|) . our 
simulations suggest that a large fraction of the gas in the 
centers of disk galaxies is ionized. 

5. CONCLUSIONS 

We have presented new MIPS images of M31 at 24, 
70, and 160 /im which have higher spatial resolution and 
sensitivity than previous infrared images. These new im- 
ages reveal the presence of two spiral arms in addition to 
a ring of star formation. The spiral arms appear to em- 
anate from a bar in the nuclear region and can be roughly 
traced out beyond the ~10 kpc ring of star formation. 
The ring of star formation forms an almost complete cir- 
cle which is offset from the nucleus and split near the 
location of M32. This ring of star formation and the dis- 
turbed spiral arms are indications that interactions may 
be affecting the morphology of M31. Dynamical models 
of M31 interacting with M32 predict morphology qual- 
itatively similar to M31's (while models with NGC 205 
instead of M32 give only small perturbations). Specifi- 
cally, both the splitting and the offset of the star forming 
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ring are seen as well as copious spiral arms. It seems that 
M31 is not an undisturbed normal spiral galaxy, but one 
which has been significantly affected by interactions. 

This paper presents a preview of the possibilities be- 
coming available to understand M31. Data covering 
the majority of this galaxy are available in the X- 
ray (XM M, Trudolyubov et al., in prep.), ultraviolet 
fOALEX. Ifhilker et al] Effll . optical (Engelbracht et 
al., in prep.), mid-infrared (Spitzer/IRAC, Barmby et 
al., in prep.), infrared (this pap er), HI (Braun et al., 
in prep.). CO ijNieten et alJl2005l) . and radio continuum 
ijBeck et al.lll989jk these will enable us to make detailed 
models on the impact that satellite galaxies have on the 
morphology of M31. 
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